MISCELLANEOUS  PAPER  HL-81-4 


ADVECTIVE  EFFECTS  IN  THE  NUMERICAL 
SIMULATION  OF  LONG  PERIOD 
LARGE  AMPLITUDE  WAVE  BEHAVIOR 

by 

H.  Lee  Butler 


Hydraulics  Laboratory 

U.  S.  Army  Engineer  Waterways  Experiment  Station 
P.  O.  Box  631,  Vicksburg,  Miss.  39180 

September  1981 
Final  Report 

I  Approved  For  Public  Mmm;  DMributioa  UnHmittd  [ 


DT!C 

■  -/Z^ECTE 


,  OCT  8  1981 


i 


.  *-0k  ■ 

'  -v-r 


lor  Assistant  Secretary  of  the  Army  (R6>D) 
Department  of  the  Army 
Washington,  D.  C.  20315 


•l*  10  7  na* 


V  v; 


i  ■ 


Destroy  this  report  when  no  longer  needed.  Do  not  return 
it  to  the  originator. 


ST,* 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated, 
by  other  authorized  documents. 


The  contents  of  this  report  are  not  to  be  used  for 
advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of 
such  commercial  products. 


t  MOMITO^INO  AOCNCV  NAME  5  ADDRCSSfl/  dllterent  from  Controlling  Office)  IS.  SECURITY  CLASS,  (el  thie  tome*) 

/  ,  _  Unclassified 


Approved  for  public  release;  distribution  unlimited 


1 17.  Distribution  statement  (at  ui.  ooouoct  Block  n,  it  dttt.m ■  >m  Knott) 


Available  from  National  Technical  Information  Service,  5285  Port  Royal  Road, 
Springfield,  Va.  22151. 

If-  KEY  WORDS  (Continue  on  rarwct  •  Ido  It  neceeeery  end  Identity  try  bloc*  number) 

Advection 

Computerized  simulation 
Mathematical  models 
Numerical  analysis 
Water  wave  periods 

10.  A  MS  TRACT  (Comf warn  Km  tmnrtm  N*  M  i.iinn,  IS  tdoottfr  ky  Hook  mnM  " 

*T  Numerical  modeling  of  water-wave  behavior  has  progressed  rapidly  in  the  last 
decade  and  is  now  generally  recognized  as  a  useful  tool  capable  of  providing 
solutions  to  many  coastal  problems.  The  U.  3.  Army,  through  various  agencies, 
has  sponsored  development  of  two-dimensional  numerical  models  for  simulating 
long  period  wave  behavior.  Many  of  these  models  have  included  the  advective 
term3  in  their  formulation.  Demonstrated  in  numerous  papers  throughout  the 
literature,  inclusion  of  advective  effects  can  lead  to  instabilities  in  the 
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PREFACE 


The  study  described  herein  was  authorized  under  the  In-House 
Laboratory  Independent  Research  Program.  All  elements  of  the 
investigation  were  conducted  at  the  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  during  the  period  November  1977  to  September 
1981  in  the  Hydraulics  Laboratory  by  Mr.  H.  L.  Butler  under  the 
direction  of  Mr.  H.  B.  Simmons,  Chief,  Hydraulics  Laboratory,  and 
Dr.  R.  W.  Whalin,  Chief,  Wave  Dynamics  Division. 

Numerical  computations  associated  with  this  work  were  performed 
on  the  CYBER  176  and  CRAY  1  computers  located  at  the  Los  Alamos 
Scientific  Laboratory  and  the  Air  Force  Weapons  Laboratory,  Kirtland 
AFB,  New  Mexico. 

Commanders  and  Directors  of  WES  during  the  course  of  the 
investigation  were  COL  John  L.  Cannon,  CE,  COL  Nelson  P.  Conover, 

CE,  and  COL  Tilford  C.  Creel,  CE.  Technical  Director  was  Mr.  F.  R. 
Brown. 
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ABSTRACT 


Numerical  modeling  of  water-wave  behavior  has  progressed  rapidly 
in  the  last  decade  and  is  now  generally  recognized  as  a  useful  tool 
capable  of  providing  solutions  to  many  coastal  problems.  The  U.  S. 

Army,  through  various  agencies,  has  sponsored  development  of  two- 
dimensional  numerical  models  for  simulating  long  period  wave  behavior. 
Many  of  these  models  have  included  the  advective  terms  in  their  formu¬ 
lation.  Demonstrated  in  numerous  papers  throughout  the  literature, 
inclusion  of  advective  effects  can  lead  to  instabilities  in  the  solution. 
This  fact  was  exemplified  in  work  of  Sloss  (1972)  in  his  attempt  to 
include  these  effects  in  storm  surge  simulation. 

The  primary  objective  of  this  investigation  is  the  assessment  of 
the  role  of  advective  terms  in  the  numerical  simulation  of  long  period, 
large  amplitude  wave  behavior.  A  major  effort  in  accomplishing  this 
objective  was  the  development  of  an  appropriate  representation  of  the 
non-linear  terms  in  the  difference  equations  of  motion.  To  this  end 
a  stabilizing  correction,  double  sweep  implicit  scheme  was  encoded 
into  an  existing  and  extensively  used  hydrodynamic  model,  WIFM.  The 
scheme  manifested  a  high  degree  of  stability  in  various  numerical  flume 
tests. 

In  previous  applications  WIFM  was  used  strictly  as  an  inland 
flooding  model.  Efforts  within  this  research  investigation  as  well  as 
concurrent  studies  permitted  extension  of  WIFM  to  include  continental 
shelf  simulation.  Comparisons  to  a  recognized  open  coast  surge  model 
were  made  and  results  indicated  WIFM  was  appropriate  for  treating  the 
global  problem  (shelf  dynamics  and  coastal  flooding)  in  a  single  grid. 

Simulations  for  the  Louisiana  coastline  under  storm  attack  by 
Hurricane  Betsy  were  made  with  and  without  including  the  non-linear 
terms'  effect.  Results  showed  that  non-linear  term  effect  on  surface 
elevation  was  minimal  but  more  noticeable  at  the  coast  near  landfall. 
Effect  on  current  patterns  was  more  pronounced.  Further  testing  is 
recommended  in  future  U.  S.  Array  surge  investigations. 
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PART  I:  INTRODUCTION 


Background 

1.  Numerical  modeling  of  water-wave  behavior  has  progressed 
rapidly  in  the  last  decade  and  is  now  generally  recognized  as  a 
useful  tool  capable  of  providing  solutions  to  many  coastal  problems. 

The  U.  S.  Army,  through  various  agencies,  has  sponsored  development 
of  two-dimensional  numerical  models  for  simulating  long-period  wave 
behavior.  Many  of  these  models  have  included  the  advective  terms 

in  their  formulation.  Demonstrated  in  numerous  papers  throughout  the 
literature,  inclusion  of  advective  effects  can  lead  to  instabilities 
in  the  solution.  This  fact  was  exemplified  in  work  of  Sloss  (1972) 
in  his  attempt  to  include  these  effects  in  storm  surge  simulation. 

2.  The  Corps  of  Engineers  has  had  to  address  the  problem  of 
providing  reliable  estimates  of  estuarine  circulation  and  coastal 
flooding  from  tides  as  well  as  large  amplitude  phenomena  such  as  storm 
surge,  tsunamis,  and  land-slide  or  explosion  generated  water  waves. 

These  simulations  are  used  to  make  sound  engineering  decisions  regarding 
the  design,  operation,  and  maintenance  of  various  coastal  projects. 

It  is  mandatory  that  the  Army  have  numerical  prediction  models  which 
can  propagate  long  period,  large  amplitude  waves  into  nearshore  and 
overland  regions.  The  success  of  base  operations,  combat  maneuvers, 
and  emergency  aid  to  civil  defense  needs,  are  linked  to  accurate 
prediction  of  catastrophic  water  waves.  The  role  of  the  advective 
terms  under  these  conditons  must  be  assessed  to  provide  a  more  complete 
understanding  of  the  associated  hydrodynamics. 

Accuracy  and  Non-Linear  Aspects  of  Difference  Schemes 

3.  The  influence  of  the  time  step  on  accuracy  of  the  difference 
scheme  is  very  important.  It  is  characterized  by  the  dimensionless 

quantity  k  3  ^  t.tV>  oro  n  -fe  flio  mav-{  ttit  im  uaira  anaoA  and  Afr  Av  ava 


where  c  is  the  maximum  wave  speed  and  At,  Ax  are 
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the  time  and  grid  step  sizes,  respectively.  When  using  explicit 
schemes,  linear  stability  investigations  require  k  <_  1.  With  implicit 
schemes,  this  restriction  usually  does  not  apply  but  accuracy  does 
diminish  with  increasing  k.  Yet  another  phenomenon  can  occur,  intro¬ 
duction  of  non-linear  instabilities  which  totally  destroy  the  solution. 
These  instabilities  have  been  shown  (Kuipers  and  Vreugdenhil,  1973) 
to  be  directly  related  to  inclusion  of  the  advective  terms  in  the 
difference  equations.  Omission  of  these  terms  makes  the  instabilities 
disappear,  but  also  makes  it  impossible  to  compute  circulation  currents 
and  horizontal  eddies.  Analysis  shows  that  flow  patterns  from  stable 
computations  are  disturbed  by  small  vortices  with  a  diameter  the  size 
of  a  grid  cell.  Vorticity  is  obviously  transferred  toward  the  smallest 
possible  scale  and  is  not  properly  removed.  The  advective  terms  are 
indispensible  for  accurate  representation  of  vorticity  transport.  The 
question  raised  in  this  report  is  whether  these  terms  are  important 
in  modeling  long  period,  large  amplitude  waves,  particularly,  in 
catastrophic  flood  situations. 


Objectives 


4.  The  primary  objective  of  this  investigation  is  the  assessment 
of  the  role  of  advective  terms  in  the  numerical  simulation  of  long 
period,  large  amplitude  wave  behavior.  An  implicit  coastal  flooding 
model,  WIFM  (Butler,  1978),  was  used  to  arrive  at  this  assessment. 

WIFM  has  been  applied  extensively  (Butler,  1981)  but  the  non-linear 
instability  problem  appeared  in  many  of  these  applications.  Thus,  a 
major  effort  in  this  study  was  to  develop  an  appropriate  representation 
of  the  non-linear  terms  in  the  difference  equations  of  motion. 

5.  Prior  to  this  study  WIFM  had  not  been  applied  to  an  open- 
coast  surge  problem.  This  investigation,  in  conjunction  with  on-going 
cost-reimbursable  studies,  provided  the  impetus  to  extend  WIFM  to 
treat  open-coast  problems.  The  objectives  were  met  by  carrying  out  the 
tasks  enumerated  below: 


a.  Develop  an  accurate  representation  of  the  nonlinear 
advective  terms  in  the  WIFM  algorithm. 

b.  Extend  WIFM  to  treat  open-coast  surge  problems  and 
compare  results  of  a  WIFM  simulation  with  a  recognized  open-coast 
model,  SSURGE  (Wanstrath,  et  al. ,  1976). 

c.  Assess  advective  effects  in  simulating  the  catastrophic 
flooding  caused  by  Hurricane  Betsy  which  struck  the  Louisiana  coast 
in  1965. 


PART  II:  COMPUTATIONAL  TECHNIQUES 


Equations  of  Motion 

6.  The  equations  of  fluid  flow  used  in  WIFM  are  derived  from 
the  classical  Navier-Stokes  equations  in  a  Cartesian  coordinate  system. 
By  assuming  vertical  accelerations  are  small  and  the  fluid  is  homo¬ 
geneous  and  integrating  the  flow  from  sea  bottom  to  water  surface,  the 
usual  two-dimensional  form  of  the  equations  of  momentum  and  continuity 
are  obtained.  A  major  advantage  of  WIFM  is  the  capability  of  applying 
a  smoothly  varying  grid  to  a  given  study  region  permitting  simulation 
of  a  complex  landscape  by  locally  increasing  grid  resolution  and/or 
aligning  coordinates  along  physical  boundaries.  For  each  direction, 
a  piecewise  reversible  transformation  which  takes  the  form 

x  =  a  +  ba^  (1) 

where  a,  b,  and  c  are  arbitrary  constants,  is  independently  used  to 
map  prototype  or  real  space  into  computational  space.  Many  stability 
problems  commonly  associated  with  variable  grid  schemes  are  eliminated 
via  the  continuity  of  the  transformation  procedure.  The  resulting 
equation  of  motion  in  a-space  can  be  written  as 


Momentum : 


.  1  .  1 

u  H - uu  +  —  vu  -  fv 

t  u1  ox  P2  a2 


+  (n  -  n  )  +  (u2  +  V2)1/2 

V1  a  ai  C2d 


/  /I  •.  2  1  /In.  ,  /I  ,2 

-e(( — )  u  +  —  ( — )  u  +  ( — )  u 

^l  alal  W1  W1  al  °1  W2  a2°2 


)  +  F 


‘2  2 


=  0 


(2) 


9 


'  r  ii-mn' 


v  +  —  uv  +  —  w  +  fu  +  (n  —  n  ) 
t  V1  *1  u2  V2  a  a2 


^  ev  ,  2  ,  2,1/2 

+  (u  +  v  ) 


1  ,1 


C2d 


(u  +  v  )  -  e(  (r~)  v  +  —  (— )  v 

V1  11  yl  yl  al  al 


+  (— )  v  +  —  (— )  v  )  +  F  =0 

^2  a2a2  ^2  ^2  a2  a2  ^2 


(3) 


Continuity: 


n„  +  —  (du)  +  — (dv) 
t  M1  a1  y2  a2 


=  R 


where 


(4) 


3x  ,  3y 

yl  =  ^  3nd  y2  =  3^ 

and  n  is  the  water-surface  elevation;  na  is  the  hydrostatic  elevation 
corresponding  to  the  atmospheric  pressure  anomaly;  u  and  v  are  the 
vertically  integrated  velocities  at  time  t  in  the  and  a2  directions, 
respectively;  d  =  n  -  h  is  the  total  water  depth;  h  is  the  still-water 
elevation;  f  is  the  Coriolis  parameter;  C  is  the  Chezy  frictional 
coefficient;  g  is  the  acceleration  of  gravity;  c  is  a  generalized  eddy 
viscosity  coefficient;  R  represents  the  rate  at  which  additional  water 
is  introduced  into  or  taken  from  the  system  (for  example,  through  rainfall 
and  evaporation) ;  and  and  Fq  are  terms  representing  external  forcing 

functions  such  as  wind  stress  in  the  and  a2  directions.  Quantities 
y^  and  u2  define  the  stretching  of  the  regular-spaced  computational  grid 
in  a-space  to  approximate  a  study  region  in  real  space.  Directions 
and  «2  correspond  to  x  and  y,  respectively.  The  single  underlined 
terms  are  referred  to  as  the  advective  inertia  terms  and  the  doubly 
underlined  terms  as  the  horizontal  diffusion  or  eddy-viscosity  terms. 
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Finite  Difference  Formulation 


7.  The  differential  equations  (Eqs.  2-4)  are  to  be  approximated 
by  difference  equations.  Various  solution  schemes,  including  implicit 
and  explicit  formulations,  could  be  used.  WIFM  permits  a  selection  of 
difference  formulations,  but  this  report  will  concentrate  on  an 
alternating  direction  technique  (ADI) . 

8.  Prior  to  the  subject  study  WIFM  employed  a  typical  ADI  scheme 
(Butler,  1978)  similar  to  Leendertse  (1970).  The  problem  with  applying 
this  procedu-''-  was  maintaining  stability  when  the  advective  inertia 
terms  were  included  in  the  solution  algorithm,  particularly  in  storm 
surge  or  long  period,  large  amplitude  wave  simulation.  Weare  (1976) 
indicated  that  the  problem  lay  in  the  differencing  techniques  used, 
namely,  in  approximating  the  advective  terms  with  expressions  not 
centered  in  time.  As  a  cure  for  this  problem  Weare  (1979)  introduced 
methods  of  developing  alternate  ADI  solution  schemes  which  did  treat 
all  terms  centered  in  time  and  space.  In  particular,  he  suggested  a 
stabilizing  correction  scheme  (SC  scheme)  employing  three  full-time 
levels.  The  details  of  this  development  can  be  found  in  a  paper  by 
Butler  (1980)  and  in  a  report  (Butler,  1981)  documenting  the  WIFM 
model . 

9.  If  the  linearized  equations  of  motion  are  written  in  matrix 
form,  one  obtains 

U  +  AU  +  BU  =  0  (5) 

t  x  y 


The  SC  scheme  for  solving  Equation  5  is 


(6) 


(1  +  Xx)  U*  =  (1  -  Xx  -  2Xy)  U*'x 
(1  +  X  )  Uk+1  =  U*  +  Xy  11*  (7) 


where 


X  «  ^  4^-  A6  and 
x  2  Ax  x 


X  =  j  ~  B6 
y  2  Ay  y 


The  quantities  Sx  and  6y  are  centered  difference  operators  and  superscript 


k  counts  time  levels.  The  starred  quantities  can  be  considered  approx¬ 
imate  values  for  corresponding  variables  at  the  (k+1)  time  level. 

10.  The  first  step  in  each  procedure  is  carried  out  by  sweeping 
the  grid  in  the  x  direction,  and  the  second  step  is  computed  by  sweeping 
in  the  y  direction.  Completing  both  sweeps  constitutes  a  full  time 
step,  advancing  the  solution  from  the  kth  time  level  to  the  (k+1)  time 
level.  The  form  of  the  difference  equations  for  the  x-sweep  is  given 
by 


Cn*  '  "k~1)  +  SJ  +Wsy  <vk"1,i)  ‘  0  (8> 

2k  <u*  -  “k'1)  +  2^  4*  <"*  +  "k*1)  ■  0  (9) 

m <v*  -  vk"1)  +  ^  5y  <"k‘1>  ■ 0  (10) 


and  the  y- sweep  by 


1  ,  k+1  .  1  .  ,  k+1,  k-1 _ 

2aT  (n  _  n*)  +  2a7  6y  (v  d  "  v  d)  ° 


(11) 


k+1  . 

u  =  u* 


1  ,  k+1  . x  g  ,  ,  k+1  k-1.  n 

2^F  (V  "  v*}  +  2A7  6y  (n  “  n  }  ° 
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1 


>1 

i 

j 


Noting  that  v*  in  Equation  10  is  only  a  function  of  previously  computed 
variables  at  the  (k-1)  time  level,  the  above  equations  can  be  simplified 
to  give 


(«*  -  nk'h  +  5X  (uk+1d  +  uk_1d)  +  -±-  6y  (vk“1d)  -  0  (14) 

BE  <“k+1  -  “k'1>  +  *  <"*  +  -'k'1>  ■  o  (15) 


1  ,  k+1  .  1  -  ,  k+1,  k-1,.  _  , , 

2At  (n  “  n  >  +  2A7  6y  (v  d  ‘  V  d)  “  0  (16) 

1  ,  k+1  k-1.  g  ,  ,  k+1  .  k-1.  „ 

IZT  (v  ■  v  J  +  6y  <n  +  n  >  "  0  (17) 

11.  Expanding  the  SC  scheme  to  the  full  equations  of  motion. 

Equations  2-4,  and  defining  the  appropriate  variables  on  each  grid  cell 
in  a  space-staggered  fashion  as  depicted  in  Figure  1, 


M 


a -FLOW/UNIT  WIDTH  IN 
X-DIRECTION  (U) 

A- FLOW/UNIT  WIDTH  IN 
Y-DIRECTION  (V) 

O-SURFACE  ELEVATION (v  ), 

WATER  DEPTH  (d), 

FRICTIONAL  COEFFICIENT  (C  OR  n) 


Figure  1.  Computational  Cell  Definition 

difference  equations  for  the  x-sweep  (along  a  grid  cell  column  parallel 
to  the  x-axis)  can  be  written  as 


1  /  .  x  ■  1  r  *  i  k+lvk  k-l-k.  . 

(n*  -  n  )  +  rr  [<$  (u  d  +  u  d  )] 


2At 


2y^Act^ 


+  Sc2  (vk“l5k)  =  R  at  <n,m) 


(18) 


1  k+1  k-lv  ,  1  k  „  ,  kx  ,  1  »k.  y  k. 

—  (u  -  u  )  +  -r - - -  u  A  (u  )  +  -x - ; -  V  d„  (u  ) 

Ip-jAc^  2axv  2y2Act2  2a2 


2At 


-f"k  +  2^^*  +  ^“1-  2^k>^ 


+  -ifr  T  uk+1  t(uk~V  +  (vk_1)2]1/2  -  e  ( — - r  6  ^  (uk) 


(C2d)k 


,  .  .2  a,  a, 
(y^Ao^)  1  1 


(y2Aa2)2  “2a2 


<“k>  +  —±-j  e  0  s  (uk) 

2y1Aa1  al  V1  /al 


+  -~2  4  0  i2  (ak)l  +  -  0  <»,»  +  i) 

2y2Aa2  222  1 


(19) 


In  these  expressions,  a  single  bar  represents  a  two-point  average  and 
a  double  bar  a  four-point  average.  The  subscripts  m  and  n  correspond 
to  spatial  locations  and  superscript  k  to  time  levels.  The  difference 
operator  6^  is  defined  as 


VZ>  '  Z«+l/2  -  Vl/2 

for  any  variable  Z. 

12.  Applying  these  equations  at  each  grid  cell  in  a  given  column 
results  in  a  system  of  linear  algebraic  equations  whose  coefficient 
matrix  is  tridiagonal.  The  y-sweep  is  formulated  in  an  analogous  manner. 
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Non-Linear  Effects 


13.  Since  the  existence  of  the  non-linear  instabilities  in 
previously  applied  difference  schemes  were  shown  to  stem  from  the 
imperfect  time-centering  of  difference  representations  of  the  non¬ 
linear  terms  (Weare,  1976),  a  fully  time-centered  scheme  was  adopted 
and  encoded  into  WIFM.  The  horizontal  diffusion  terms  also  were 
included  in  the  difference  scheme  to  contribute  to  numerical  stability 
(Vreugdenhil,  1973).  These  terms,  strictly  speaking,  should  be 
included  in  the  momentum  equations.  Vreugdenhil  demonstrates  that 
such  terms  actually  are  representations  of  the  effective  stress  terms 
usually  neglected.  As  discussed  in  following  paragraphs  tests  were 
made  to  determine  effects  of  each  term. 


Extension  To  Open-Coast  Treatment 


14.  In  previous  surge  applications  (Butler  and  Wanstrath,  1976); 
Butler,  1978)  WIFM  was  used  strictly  as  an  inland  flooding  model. 
Seaward  boundary  conditions  were  obtained  from  open-coast  surge  models 
by  specifying  surface  elevation  at  the  boundaries.  To  extend  WIFM's 
capabilities  to  include  shelf  simulation  only  treatments  of  two 
additional  boundary  conditions  were  added.  These  were  a  uniform  flux 
boundary  condition  and  an  inverted  barometer  effect  condition. 
Capability  of  using  the  inverted  barometer  effect  throughout  the  grid 
for  initial  conditions  also  was  added. 


PART  III:  TEST  APPLICATIONS 


Flume  Tests 


15.  To  test  the  new  WIFM  algorithm,  including  representation  of 
both  advective  and  horizontal  diffusion  terms,  two  simple  flume  tests 
were  devised.  Both  tests  used  a  constant  depth  flume  of  10  ft*.  One 
test  was  for  an  expanding  flow  and  the  second  for  flow  around  a  break¬ 
water.  Various  cell  sizes  were  employed  to  test  spatial  resolution 
effects.  Two  types  of  assumptions  were  made  for  approximating  the 
non-linear  terms  at  closed  boundaries:  linerization  at  the  boundary 
(Al)  or  one-sided  difference  quotients  to  replace  these  terms  (A2) . 
Linerization  introduces  additional  boundary  conditions  approximating 
free  slip  which  may  influence  secondary  phenomena. 

16.  A  constant  value  for  the  eddy-viscosity  coefficient,  e,  was 
assumed  over  the  spatial  regime.  Various  values  for  e  were  tested. 
Increasing  £  beyond  physically  justifiable  limits  seriously  affected 
both  surface  elevation,  but  more  particularly,  the  flow  pattern. 

Figures  2  and  3  depict  circulation  patterns  for  an  expanding  flow 

for  480  time  steps  and  assumptions  Al  and  A2,  respectively.  A  uniform 
flow  was  imposed  at  the  left-hand  boundary.  Figures  4  and  5  depict 
full  development  after  960  time  steps.  Boundary  effects  for  assumption 
Al  are  obvious  when  a  steady  state  condition  was  reached  (Figure  4) . 

17.  A  barrier  was  inserted  in  the  flume  and  a  uniform  flow  was 
again  imposed  at  the  left.  Figures  6-12  display  sample  tests  for  varying 
spatial  resolution.  Eddys  are  developed  behind  the  breakwater  as 
expected.  When  the  resolution  was  decreased  to  only  4  cells  across  the 
tank  no  eddy  developed.  Figures  11  and  12  indicate  that  even  a  short 
barrier  extension  of  two  cell  lengths  will  produce  an  eddy.  Assumption 
Al  and  A2  produced  almost  identical  results  except  for  slight  differences 
at  the  tip  of  the  breakwater.  No  definitive  conclusion  can  be  made  on 
which  assumption  is  better,  but  assumption  A2  is  a  more  accurate  model  of 
the  physics  involved  and  thus  should  be  used. 

*  Conversion  factor  for  feet  to  metres:  multiply  by  0.3048 
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Figure  3.  Expanding  flow  for  assumption  A2  after  480  time  steps 
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Figure  9.  Enlargement  of  Figure  8  in  area  of  barrier 
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Surge  Applications 

18.  Since  WIFM  capabilities  were  extended  to  include  modeling  of 
the  continental  shelf  a  comparison  study  between  it  and  a  recognized 
open-coast  surge  model  was  performed.  The  basis  for  the  comparison  was 
the  problem  of  hindcasting  water  elevations  caused  by  Hurricane  Betsy 
(September  1965).  SSURGE,  an  open-coast  model  by  Wanstrath  (1978), 
was  previously  applied  to  this  problem  in  a  study  of  coastal  flooding 
for  Louisiana.  The  windfield  model  used  was  a  slight  modification 

of  an  analytical  model  formulated  by  Jelesnianski  (1965).  All  windfield 
parameters  from  the  previous  study  were  retained  for  use  in  the 
present  comparison  effort.  Land  effects  on  the  windfield  were  not 
simulated.  Figures  13  and  14  depict  the  Wanstrath  and  WIFM  grids  used. 

Figures  15  and  16  display  hydrograph  results  at  two  representative 
gage  loccations.  SSURGE  does  not  account  for  flooding  but  does 
permit  ponding  areas  through  a  finite  height  barrier  coast.  The  grid 
extended  only  to  the  300-ft  contour.  The  WIFM  grid  extends  to  depth 
contours  greater  thean  2,000  ft  and  treats  the  two-dimensional  overland 
flooding  and  hence  should  be  able  to  portray  the  surge  more  accurately. 

The  mean  absolute  error  in  predicting  open-coast  observed  peak  watersurface 
elevations  was  essentially  the  same  for  both  models  (approximately  0.7  ft). 

19.  The  WIFM  grid  shown  in  Figure  14  was  devised  especially  for 
the  comparison  study.  A  grid,  extending  over  a  larger  coastline  reach 
and  containing  higher  resolution  (shown  in  Figure  17) ,  was  constructed 

for  use  in  an  ongoing  study  for  the  U.  S.  Army  Engineer  District,  New  Orleans. 
Simulation  of  Betsy  was  repeated  but  this  time  a  version  of  the  Standard 
Project  Hurricane  (SPH)  parametric  windfield  model  was  used,  one  which 
introduces  deformations  to  the  windfield  to  account  for  land  effects. 

This  model  of  the  coast  was  run  with  and  without  the  advective  and 
horizontal  diffusion  terms.  Effects  (on  surface  elevations)  of  including 
non-linear  terms  in  the  computations  were  minimal  throughout  the  model 
region.  The  only  noticeable  effect  was  in  the  coastal  area  near  landfall 
(Grand  Isle,  LA)  where  surface  elevations  were  0.3  to  0.6  ft  higher. 
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Figure  13.  SSURGE  computational  grid 
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Figure  15.  Comparison  of  WIFM  and  S SURGE  hydrographs  with  observed 
data  at  Biloxi,  MS 


V 


*  V.' 


Figure  17.  WIFM  open-coast  computational  grid  for  the  Lake  Pontchar train 
Hurricane  Barrier  Study 


This  increase  in  water  level  was  a  result  of  stronger  non-linear  term 
influence  on  circulation.  Surface  elevations  generally  increase  near 
the  tip  of  the  Mississippi  River  Delta  (to  the  right  of  landfall) . 

Figure  18  depicts  elevation  and  velocity  comparisons  for  a  gage  located 
at  the  tip  of  the  Delta.  Changes  in  velocity  mangitude  and  direction 
at  this  gage  are  typical  of  near-coast  circulation  differences. 

20.  Discharge  ranges  were  established  throughout  the  model  area. 

A  2%  reduction  of  transport  through  the  major  artery  to  Lake  Pontchartrain 
was  noted.  Most  other  range  comparisons  showed  the  same  minimal  effect. 
The  comparison  for  a  range  aligned  north-south  across  the  eastern  end 
of  the  lake  exhibited  the  largest  difference  in  discharge  (25%  more  in 
the  ebb  direction  with  the  non-linear  terms  included) .  Since  the  range 
is  long  (about  8  miles)  and  circulation  is  unrestricted  the  effect  of 
the  non-linear  terms  on  circulation  is  expected  to  be  more  significant. 
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PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


21.  An  accurate  method  for  representing  the  non-linear  terms  in 
the  equations  of  motion  was  developed.  No  numerical  stability  problems 
were  encountered  while  running  the  various  tests  described  in  this 
report.  Additional  tests  for  varied  difference  approximations  of 
these  terms  near  land-water  boundaries  are  recommended  for  future 
research.  In  particular,  tests  of  free-slip  and  no-slip  representations 
are  suggested.  An  attempt  was  made  in  this  study  to  look  at  this 
problem  but  no  definitive  conclusions  were  reached.  No-slip  represen¬ 
tations  tested  indicated  much  too  severe  of  an  effect  on  the  current 
magnitudes  in  the  cell  adjacent  to  the  boundary.  WIFM  was  extended  to 
simulate  open-coast  surge  on  the  continental  shelf  and  a  successful 
comparison  for  Hurricane  Betsy  to  a  recognized  open-coast  model,  SSURGE, 
was  made.  This  development  permits  shelf  and  overland  surge  simulation 
on  the  same  grid  substantially  reducing  computational  costs  in  future 
applications. 

22.  Research  funding  did  not  permit  testing  for  advective  effects 
in  surge  computations  over  a  wide  variety  of  conditions,  both  meteoro¬ 
logical  and  varied  land  configurations.  Catastrophic  flooding  from  Betsy 
was  simulated  with  and  without  the  non-linear  terms.  Results  from  these 
runs  showed  that  minimal  effects  on  water-surface  elevations  and 
current  magnitudes  were  caused  by  including  the  non-linear  terms.  The 
largest  effect  was  on  the  peak  surge  at  the  coast  near  landfall  (increases 
up  to  0.6  ft)  and  circulation  near  the  coast.  The  total  discharge  through 
the  major  lake  entrance.  The  Rlgolets,  was  decreased  by  a  little  over 

2%.  Based  on  this  study  the  non-linear  terms  could  be  neglected  if 
interests  were  solely  water  surface  elevations.  On  the  other  hand,  if 
current  patterns  were  required  (for  example,  to  move  sediment,  pollutants, 
etc.)  the  results  indicate  that  these  terms  should  be  included.  Further 
testing  is  recommended  in  future  U.  S.  Army  surge  investigations. 
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